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1975 Annual Report on the Investigation 
of Critical Burning of Fuel Droplets 


Summary 


This report discusses activities under NASA Contract NCR 39-009- 
077 for the period January 1, 1975 to December 31, 1975. During this 
period the work was divided into two phases: (1) Combustion response 
of monopropellant droplets and sprays and (2) High pressure combus- 
tion of bipropellant sprays. The accomplishments of each phase of 
the study for the present report period may be summarized as follows: 

1. Combustion response of monopropellant droplets and sprays . 
Earlier results for the response of burning liquid mono- 
propellant strands were extended to consider monopropellant 
droplets and sprays. The theory employs perturbation 
analysis with calculations through first order to yield the 
linear response. Physical properties were chosen to simu- 
late hydrazine decomposition due to the availability of 
data for this material to test the model. The zero-order or 
steady state calculations gave predictions in good agreement 
with earlier measurements of liquid temperatures and burning 
rates. A simple, approximate formula was developed which 
estimates actual droplet burning rates with a maximum error 
of lOZ; this formula is readily Integrable and convenient 
for use in the analysis of spray combustion. 

The response of individual droplets was examined. Small 
droplets gave low or negative response at all frequencies. 
Large drops provide a response near unity at low frequencies, 
with the response declining at frequencies somewhat in excess 
of the characteristic frequency of the liquid phase thermal 
wave. If temperature gradients are present in the liquid 
phase, this decline can be preceeded by a response peak in 
excess of unity. A second response peak is also present 
for large droplets which is associated with gas phase 
transient effects, however, the frequency range of this 
peak is generally higher than the frequ'sncy range normally 
associated with combustion instability. 

The response of a monodisperse injected spray is less than 
the response of the initial droplet size, due to the re- 
duced response of small droplets near the end of their 
lifetime. The response characteristics of a polydisperse 
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Injected spray depends strongly on the size distribution and 
several cases were considered. Notably, even a very low per- 
centage o£ large droplets can yield appreciable response, due 
to their relatively long lifetimes and high response. 

Work In this area of the Investigation has been concluded. 

High pressure spray combustion . An apparatus has been de- 
signed and fabricated to allow observation of fuel spray 
combustion at pressure levels where thermodynamic critical 
phenomena are Important. The apparatus provides for spray 
shadowgraphs, dark field flame photographs and gas tempera- 
tures In the region downstream of the spray region. 

Modeling of this process Is based on an Integral-entrainment 
model which allows for density variation In the two-phase 
flow and high pressure thermodynamic characteristics, but 
neglects slip between the phases. This type of model is most 
accurate at high pressures and Initial theoretical and experl 
mental comparisons are considering this regime. 

The assembly of the apparatus is nearing completion and test- 
ing will begin early in the next report period. Initial 
development of the theory has also been completed. 
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1975 Annual Report on the Investigation 
of Critical Burning of Fuel Droplets 


1. Introduction 


The objective of this Investigation was to continue earlier work 
on the steady and unsteady combustion of liquid fuel droplets under 
rocket engine conditions. Emphasis has been placed on consideration 
of combustion at elevated pressures and temperatures, representative 
of typical combustion chamber environments. The results of this study 
have technical application to the design of liquid fueled rocket engines 
and the determination of the combustion Instability characteristics of 
these engines. 

This report gives a summary of progress on the Investigation for 
the period January 1, 1975 to December 31, 1975. During this report 
period, the work was divided Into two phases, as follows: 

1. Combustion response of monopropellant droplets and sprays. 

This portion of the study Is concerned with extending the 
analysis of the "open loop" response of burning liquid mono- 
propellant strands (1,2), to the more practical case of drop- 
lets and sprays. The theoretical approach employs a perturba- 
tion analysis allowing for both liquid and gas phase transient 
effects. For conditions of technological Interest, the gas 
phase was found to be largely quasisteady and present calcula- 
tions have emphasized this condition. The analysis of spray 
response was limited to the quasisteady gas phase approxima- 
tion. The procedure of the study Involved considering the 
steady combustion and the response of Individual droplets 
first, follows I by analysis of the steady combustion and re- 
sponse of spr/^ys. Work In this area was concluded during the 
report period. 

2. High pressure spray combu tlon. An investigation was begun 
on the characteristics of burning fuel sprays at elevated 
pressures where critical phenomena are Important. The 
modeling of this process is based upon an Integral-entrainment 
model, with provision for Including the high pressure thermo- 
dynamic characteristics of the liquid gasification process. 
During this report, an apparatus to test the theory was de- 
signed and fabricated, and Initial development of the theory 
was begun. 

Reports and papers issued during the period of this report are 
as follows: 

G. M. Faeth and S.Chanin, "1974 Annual Report on the Investi- 
gation of Critical Burning of Fuel Droplets," Contract NCR 39- 
009-077, The Pennsylvania State University, University Park, 

PA, February 1975. 
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C. B. Allison and G. M. Fssth, "Open-Loop Response of a 
Burning Liquid Monopropellanc • " AlAA J 13,10,1287-1294 (1975). 

6. S. Canada and G. M. Faech, "Combustion of Liquid Fuels In 
a Flowing Combustion Gas Environment," Fifteenth Symposium 
(International) on Combustion, pp. 419-428, The Combustion 
Institute, Pittsburgh (1975). 
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2. Coabuitlon Rtiponst of Monoproptllant Droplets and Sprays . 
2.1 Introduction 


Th« objective of this Investigation Is to study the burning 
rate response of monopropellant droplets to Imposed pressure 
oscillations, followed by consideration of the total response of 
a spray. Knowledge of the frequency ranges where the combustion 
response Is high, allows the designer either to adjust combustion 
geometry to avoid having characteristic chamber frequencies In 
this region, or provide damping tuned to this frequency range. 

The present work Is a continuation of earlier efforts which 
created the response of a burning monopropellant strand to Imposed 
pressure oscillations (1,2). Using these findings as a basis, the 
.iiore practical case of droplet combustion was Investigated. By 
utilizing a model generally substantiated by strand combustion re- 
sults, the case of droplet combustion could be approached with some 
confidence. 

In the following, only a brief summary of the oscillatory 
combustion of monopropellant droplet work Is given. A complete 
discussion of this Investigation may be found In Reference A. 

2.2 Theory 

The theoretical model under consideration Involves a droplet 
that Is burning In the absence of convection, yielding a spherical- 
ly symmetric flow field. While convection Is Important under act- 
ual combustion chamber conditions, there Is evidence to Indicate 
that the region where the droplet has a low (or zero) relative 
velocity with respect to the ambient gas, Is a very critical zone 
In determining combustion Instability characteristics (5). For 
simple evaporation without decomposition, the response due to 
velocity effects has been found to be quite low (6). In addition, 
the decomposition process near the droplet surface reduces the 
Influence of convection over a fairly vide range of conditions 
(7,8). This occurs since convection only affects the process 
when the outer edge of the flow field Interacts with the reaction 
zone, a situation that is limited to very weak reactivity or very 
high Reynolds numbers. Therefore, neglecting convection for mono- 
propellant droplets puts fewer limitations on the practicality of 
the calculations, than would be the case for bipropellants. 

Similar to earlier response studies (1,2,6), the ambient 
pressure is assumed to be oscillating with a wavelength that is 
long In comparison to the dimensions of the combustion field of 
the droplet. The period of oscillation, however, is assumed to 


4 


bft short In conpsrlson to th« total droplet llfotlna to that 
largo changoa In tho position of tha surfaca during a prassura 
oscillation ara not considared. This assumption allows tha 
analysis to procaad vhila only examining oscillatory solutions. 

Tha low frequency regime where this assumption breaks down is 
associated with the range of frequencies where both liquid and 
gas phases are quasi-steady. 

Since the combustion rate of monopropellant droplets is much 
higher than bipropellants, the usual assumption of a constant 
liquid phase temperature (in the zero order) is less valid than 
for bipropellants. Therefore, the presence of mean liquid phase 
temperature gradients is considered in the analysis. Examination 
of constant mean temperatures is also made, in order to Include 
the conventional steady combustion model, by equating the bulk 
liquid temperature to the wet bulb temperature at a given pressure. 

A sketch of the theoretical model is shown in Fig. 1. The 
process is examined at an Instant of time when the droplet radius 
is r* (quantities with an * are dimensional). Formally, this 

radius is taken to be fixed so that the mass flux of fuel is time 
varying. This actually corresponds to porous sphere combustion; 
however, the two cases are equivalent as long as the density in 
the liquid phase is large compared to that of the gas phase; and 
the period of oscillation is small in comparison to the lifetime 
of the d) >let (9). When this is true, the motion of the surface 
with rest :t to the mean surf ce position is negligible and can 
be Ignored. Exceptions to th.s assumption arise near the critical 
point and the present analysis is not valid in this r*<riime. The 
response portion of the analysis is invalid at frequencies having 
an oscillation period comparable to the lifetime of the droplet. 

For generality, gas phase transient effects are included in 
the model. The effect of variable properties Is also Included, 
so that the model is equivalent to the earlier strand combustion 
analysis (1,2). Hydrazine was the fuel studied, since earlier 
information was available for both oscillatory (1,2) and steady 
droplet burning (10) characteristics. 

The remaining assumptions of the analysis are similar to 
those of Refs 1 and 2. They are as follows: 

1. The process is spherically symmetric with a Mach number 
much less chan unity and negligible body forces. Inertial and 
viscous terms in the momentum equation are neglected. 

2. The reaction process is premixed and laminar. A one-step, 
irreversible chemical reaction cakes place in the gas phase and 
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«ny cln« lags assoclaCtd wlch the cnemlcal reaction Itself are 
negligible, l.c. the chemical reaction Is locaTly quasi-steady and 
obeys an Arrhenius equation. Chemical reaction Is neglected In 
the liquid phase. 

3. Radiation heat transfer is neglected. 

4. The gas phase Is taken to be an ideal gas and the Levis 
number Is assumed to be unity. 

5. All gas diffusion coefficients are equal, all molecular 
weights are equal, all gas phase specific heats are equal and con- 
stant, the gas phase thermal conductivity Is Independent of com- 
position and varies linearly v'lth temperature, and the liquid is 
composed of a single chemical species having constant properties. 

6. The combustion products are assumed to be insoluble In 
the liquid phase and the gas phase fuel mass fraction at the 

1 it' lid surface Is given by the Clauslus-Clapeyron equation. As 
in Che case of the gas phase reaction, the equilibrium at the 
curface is assumed to occur rapidly in comparison to other 
transient effects In the system. 

7. The wavelength of any periodic pressure disturbance is 
assumed to be long compared with the diameter of the zone in- 
volving active combustion. Consideration of the momentum equa- 
tion, along with Assumption 1, then Implies that pressure is only 
a function of time. 

A discussion of the applicability of these assumptions is 
provided in Refs. 1 and 2. 

Utilizing the above assumptions, the conservation equations 
and boundary conditions are written in spherical coordinates. 

After nondimensionalizing, a perturbation solution is used with 
the amplitude of the pressure oscillation being the perturbation 
parameter, C ('t small quantity). The resulting equations are 
separated into like powers of z. The zero order problem de- 
scribes the steady state situation, and the first order problem 
describes the effect of oscillating pressure. Reference 4 
should be consulted for the details of the analysis. 

2. 3 Steady State Results 


The physical property values and the chemical kinetic param- 
eters used in the calculations are summarized in Ref. 2. 

A second order reaction correlates the strand burning rate 
of hydrazine at pressures above 1 atmosphere (1,2). For steady 
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state conditions, the theoretical burning rate is relatively 
insensitive to variations in the activation energy. However, for 
oscillatory combustion, the activation energy influences both the 
amplitude and phase of the combustion response of strands. Alli- 
son and Faeth (2) found that the value of E* (activation energy) 
gave good agreement between theoretical and experimental response 
determinations over their test rangt.. 

In order to check the accuracy of the zero-order model, pre- 
dicted and measured liquid surface temperatures were compared as 
a function of pressure. These results are shown in Fig. 2. The 
agreement is good throughout the pressure range of the measure- 
ments (.51-19.8 atm). 

The test conditions shown in Fig. 2 involved a liquid temp- 
erature of 298®K far from the liquid surface. This is representa- 
tive of a non-wet bulb condition, with the liquid temperature 
increasing as the surface is approached. A theoretical calcula- 
tion for the wet bulb state, where the liquid temperature is con- 
stant, is also presented in Fig. 2. For the wet bulb case, no 
energy is required to raise the liquid to the vaporization temper- 
ature, resulting in a higher surface temperature than the non-wet 
bulb case, for a given pressure. 

Steady state calculations were completed over a range of A 
(dimensionless parameter proportional to droplet diameter squared) 
and pressure for both non-wet bulb (with a centerline droplet 
temperature of 298°K) and wet bulb conditions. 

The mass burning rate of a liquid droplet depends on the 
size of the droplet (at a given pressure) and approaches two 
asymptotic limits: the large drop and small drop limit burning 
rates. The small drop limit represents a droplet vaporizing 
without reaction, and the large drop limit approximates a one- 
dimensional liquid strand with the flame very close to the liquid 
surface. 

Figure 3 shows the dimensionless steady state mass burning 

rate per unit solid angle, m , as a function of A at 1 atmosphere 

“0 

pressure for non-wet bulb conditions. The complete steady state 

solution clearly approaches the large and small drop limits. In 

order to plot m as a function of A at the large drop limit, 

^ 0 

values of the burning rate eigenvalue were obtained from the strand 
combustion calculations, (1,2) and converted to the present 
notation. The small drop limit mass burning rate, m , is 

SDL 
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Figure 2 Theoretical and Experimental Liquid Surface 
Temperature as a Function of Pressure 
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SMALL DROP LIMIT 


Figure 3 Tlieoretical and Experimental Steady State Droplet Gasification Rates at 
Approximately ] atm. Pressure (Non-Wet Bulb Conditions) 
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calculated directly from the theory for A-0. 

An approximation to the complete steady state solution Is 
shown in Fig. 3. The approximation was obtained by simply add- 
ing the small and large drop mass burning rates at each value of 
A, as follows: 

m ■ m + m 
o o o 

SDL LDL 

This approximation yields a maximum error of approximately 12 % 
over the entire range of A, and it will be used later to evaluate 
spray combustion characteristics since it can be readily integrat- 
ed to determine droplet lifetimes. 

Data on hydrazine combustion was available from experimental 
work done by Allison and Faeth (10) . This study Involved droplet 
burning in a combustion gas under decomposition conditions, for 
various droplet sizes, at atmospheric pressure. The data was 
limited to the ambient temperature range 1660-2530°K; therefore, 
it was necessary to extrapolate the measurements to the 13A5®K 
ambient temperature which correspond to adiabatic decomposition 
for an initial fuel temperature of 298K. The agreement between 
the predicted and measured burning rates is seen to be reasonably 
good, lending confidence to the theoretical model. 

Burning rate results for wet bulb conditions are similar to 
the non-wet bulb conditions; the major difference involves a slight 
increase in the burning rate throughout the entire range of A. 
Theoretical mass burning rates at pressures of 10 and 100 atmos- 
pheres also gave results very similar to Fig. 3 for both wet and 
non-wet bulb conditions. In all cases the transition region 
between the }.arge and small droplet limits falls approximately in 
the range 10‘*<A<109. 

In order to employ the approximate burning rate expression, 
it is necessary to have values for m and the strand burning 

°SDL 

rate eigenvalue, A^. These quantities, along with surface temp- 
erature predictions, are summarized in Table 1 over a range of 
pressures, for both wet bulb and non-wet bulb conditions. 

With increasing pressure, the liquid surface temperature, 

T *, increases in all cases. As noted earlier, the wet bulb 
os 

surface temperature slightly exceeds the non-wet bulb surface 
temperature at each pressure. 
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Table 1 

Summary of Steady State Calculations 


*(N/m^) 

OS 

«o 

m 

°SDL 

A 

s 

Non Wet 

Bulb Conditions 



10^ 

374 

1345 

.544 

1.432X10^ 

io‘ 

450 

1345 

.509 

1.432X10^ 

10^ 

562 

1345 

.454 

1.432X10^ 

Wet Bulb Conditions 



10^ 

376 

1417 

.603 

6.205X10^ 

10^ 

457 

1498 

.629 

2.660X10^ 

10^ 

582 

1623 

.668 

8.162X10^ 


center temperature 298K. 
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While A Is relatively constant with pressure at non-wet 
s 

bulb conditions, A decreases with Increasing pressure for the 
8 

wet bulb case, which significantly Increases the mass burning 
rates for the larger drops at a given value of A. This effect 
Is due to Increased reaction rates for the higher flame temper- 
atures of wet-bulb combustion, at elevated pressures. 

At non-wet bulb conditions, m decreases with Increasing 

°SDL 

pressure due to the reduction In the temperature difference 
between the ambient gas and the liquid surface. For the wet 
bulb case, the Increased ambient temperature compensates for 
this effect and there Is a slight Increase In m^' with Increas- 
ing pressure. SDL 

Using the properties from Ref. 4, the droplet size range 
corresponding to the transition region (taken to be 10^<A<10^) 

Is shown In Table 2. As the pressure Increases, the droplet 
size for the onset of the large drop limit Is reduced. Notably, 
the size range of technological Importance for actual combustors 
falls largely In the transition region. At higher pressures, 
however, a greater percentage of the droplets present In a 
spray can be represented by large drop limit results. 

2.4 Droplet Response 

In order to analyze response, the burning rate response 
function, P^, Is utilized (5). In the present notation this 

quantity Is 

P - Re {m. /m } 
r 10 

where Re denotes the real part; l.e. that portion of the burning 
rate oscillation, that is in-phase with the pressure oscilla- 
tion. For instability, P^ must be positive and of order unity 

(the exact value depends upon the degree of damping present) at 
a point in the combustion chamber where the pressure is varying. 
It must be understood that a droplet with P^ ^ 1 will not in it- 
self cause instability, but rather the sum of the responses of 
all the droplets within the combustion chamber must be greater 
than the available damping, in order for the combustion oscilla- 
tion to grow. 

At low frequencies for the small droplet or evaporative 
limit, the response function is negative, with a relatively small 
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Table 2 

Droplet Size Range of Transition Region 


* +5 

P X 10 

(N/m^) 

* 

d 

" d * (U) 

0 

®mln 

s 

max 

1 


25 - 

• 7500 

10 


2.5 

- 750 

100 


.25 

- 75 


a A 9 

Corresponding to 10 <A<10 
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magnitude and with Increasing pressure, the response function 
becomes more negative (4). This Is In contrast to results for 
bipropellant combustion where the response function generally 
approaches zero or a positive limit for small u). 

Moving from the evaporative limit (small drop limit) to the 
large drop limit, much greater response Is observed for mono- 
propellant droplets. Figure 4 Illustrates the non-wet bulb 
response as a function of frequency u (dimensionless frequency 
for strand combustion from Refs. 1 anl 2 for the large droplet 
limit at a mean pressure of 10 atmosphere. Figure 5 shows 
similar results at wet bulb conditions. In addition to the 
quasi-steady gas phase approximation, results for the completely 
unsteady gas phase are also shown. The non-wet bulb results are 
replotted from the calculations of Ref. 1, while the wet bulb 
results were calculated during the present Investigation. 

In all cases, the response approaches unity at low fre- 
quencies, where both liquid and gas phases are quasi-steady. 

This follows from the fact that the burning rate Is proportional 
to pressure, under steady conditions, for a second-order reaction 
( 1 ). 


With Increasing frequency, at non-wet bulb conditions, a 
peak Is observed In the response plot at frequencies near the 
characteristic frequency of the liquid phase. Beyond this peak, 
the quasi-steady gas phase solution gives a continuously de- 
clining response. The peak Is absent In all bulb cases, with 
the response showing a notlcable decline at the liquid phase 
characteristic frequency. This Interesting contrast between 
the two cases will be discussed later In more detail. 

At higher frequencies, the analysis allowing for transient 

gas phase effects begins to diverge from the quasi-steady gas 

phase solution. Substantial differences between the^two models 

are encountered for values of u on the order of lO”"^, which 

s 

represents frequencies near the characteristic frequency of the 
gas phase. At dimensionless frequencies on the order of unity, 
a second response peak is observed, in all cases, for the com- 
plete unsteady gas phase solution. This peak is clearly associ- 
ated with gas phase transient effects and yields maximum values 
of on the order of unity. With increasing pressure, the 

frequencies of the two response peaks approach one another due 
to the increase in gas density. 

Table 3 lists the frequencies where lar;^e droplets have 
response peaks. In considering this table, it should be recalled 
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Figure A Response as a Function of Frequency at the Large Droplet Limit at 
Approximately 10 atm. Pressure (Non-Wet Bulb Conditions). 
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Figure 5 Response as a Function of Frequency at the Large Droplet Llnit 
at Approximately 10 atm. Pressure (Wet Bulb Conditions). 
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Table 3 


Frequencies of the Unsteady Gas and Liquid Phase Response Peaks 


Mean 

Pressure 

(N/m^) 

Liquid Transient Peak (Hz) 

Gas Transient Peak (Hz) 

lo’ 


(Non-Wet Bulb Condition) 
0.U4 

528 

io‘ 


17 

11,000 

lo' 


1100* 

110,000 

10^ 


(Wet Bulb Condition) 

€70 

io‘ 


— 

6,700 

10^ 


— 

67,000 


^Corresponds with first gas phase transient analysis peak 
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chat acoustic Instability In rocket engine combustion chambers Is 
generally associated with Che frequency range 300-30,000 Hz (5). 

Ac pressures below 10 atmospheres, response peaks due to gas phase 
transient effects largely fall within this range. At pressures 
above 10 atmospheres, response peaks due to gas phase cranr^unc 
effects decline In Importance since they are associated wlr! viry 
high frequencies. In this pressure range, liquid phase phenomena 
are more significant, providing response peaks In the critical 
frequency range. Since the pressure range usually encountered In 
rocket engines is above 10 atmospheres, gas phase transient 
effects do not appear to be a major factor In causing Instability 
for hydrazine fueled engines. Based on this finding, further 
analysis In Che transition region, between the large and small 
droplet regimes, was limited to the case of quasi-steady gas phase 
effects. 

Figure 6 shows Che quasi-steady gas phase non-wet bulb 
response as a function of frequency, u, (which Is proportional 
to u)^) at various drop sizes, ranging from Che small drop to the 

large drop limit. As A, or drop size. Increases, the response 
Increases, caused by the combustion zone moving closer to the 
liquid surface. At an A of approximately 10^, a peak In the 
response curve develops. Increasing In magnitude and frequency 
as the drop size Increases, until the large drop approximates 
the one-dimensional strand near the value of A"10^^. As the 
pressure Increases, the peaks Increase In magnitude and move to 
higher frequencies. 

The quasi-steady gas phase wet bulb responses are plotted in 
Fig. 7, as a function of frequency for the complete range of drop 
sizes at 10 atmospheres pressure. As in the non-wet bulb case, 
the response increases with Increasing A; however, no response 
peaks are observed. Droplet combustion response Is well approxi- 
mated by the one-dlmenslonal stra'.d at values of A greater than 


The only two types of monopropellant droplets which have 
been considered, are the vet bulb droplet with a zero mean 
temperature gradient (constant mean temperature) and the non-wet 
bulb droplet with a steep mean temperature gradient. For the 
parameters used, a response peak Is only observed In the non-wet 
bulb case. Under realistic conditions, the droplets will have 
temperature gradients ranging between the wet and non-vet bulb 
limits. 

From the present study it can be concluded that a response 
peak will only occur when the mean temperature gradient near the 
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Figure 6 Quasi-Steady Gas Phase Response as a Function of Frequency and / 
at Approximately 10 atm. Pressure (Wet Bulb Conditions). 
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droplet surface has sufficient steepness so that enough energy 
will flow from the surface towards the center of the drop, caus- 
ing the perturbation surface temperature to lead the pressure. 

When the mean temperature gradient Is not steep enough, the per- 
turbation surface temperature will never lead the pressure. Also, 
as the droplet size decreases while It Is reacting, the liquid 
temperature gradient will tend to approach the wet bulb limit; 
thus, eliminating the response peak. 

2.5 Spray Response 

A combustion chamber will become unstable when the total 
response of all the droplets In the spray exceeds the available 
damping. Using the droplet response results, an analysis of 
spray response was undertaken. A one-dlmenslonal spray combustion 
model was used. In which the pressure field Is Identical for all 
parts of the spray. 

Droolst lifetimes were calculated using the steady state 
mass burning rate approximation. Table 4 lists the droplet life- 
times as a function of pressure and Initial drop radius. 

The wet bulb droplets consistently have shorter lifetimes 
because of their higher mass burning rates. For all cases, at a 
given Initial radius, the lifetime decreases with Increasing 
pressure. This is caused by the fact that a given radius cor- 
responds to a larger value of A at higher pressures, yielding 
a greater mass burning rate. At a constant pressure, the life- 
time Increases with drop size, which to be expected. 

The response portion of the present analysis is not valid 
at frequencies having an oscillation period comparable to the 
lifetime of the droplet. This provides a minimum droplet size 
that can be considered by the present analysis. The probable 
oscillation periodr. of instability range from 2 x 10”^ to 

3 X 10 seconds/cycle (corresponding to the 500-30,000 Hz range 
indicated in Ref. 5). At the oscillation period limit of 

_3 

2 X 10 seconds/cycle, the minimum drop radius is approximately 
13 u at the non-wet bulb condition, and approximately 95 U for 

the wet bulb case. At tne other limit of 3 x 10 ^ secor.ds/cycle 
the minimum droplet radius is approximately 2 y at all pressures 
and both wet and non-wet bulb conditions. This indicates that at 
low frequencies, appreciable response may be generated py droplet 
lifetime characteristics, a mechanism that has not been considered 
in the present investigation. 


Table 4 


Dimensional Droplet Lifetimes as a Function 
of Initial Radius and Pressure 


P„*(N/m^) r *(U) “ 10 100 1000 

o so 


(seconds) 


Non-Wet Bulb 




10^ 

0.00158 

0.128 

4.79 

10^ 

0.00135 

0.0489 

0.728 

10^ 

0.000507 

0.00734 

0.0788 

Wet Bulb 

(seconds) 


10^ 

0.00141 

0.109 

3.61 

O 

1 O' 

0.000955 

0.0272 

0.361 

10 

0.00179 

0.0020 

0.0229 
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If a monodlsperse spray Is Injected Into the combustion 
chamber, the total response, P , Is the weighted sum of all drop 
sizes existing In the chamber.^ Figure 8 Illustrates the non-wet 
bulb total response of a monodlsperse spray as a function of 
frequency, u , and (dimensionless parameters based on the 

Injected size) for a pressure of 10 atmospheres. The solid lines 
correspond to the entire spray while the dashed lines correspond 
to the situation when It Is assumed that all drop sizes are 
approximated by a constant drop size corresponding to A^. For 

these two cases, the peaks occur at approximately the same 
frequency, while the magnitude of the peak Is slightly larger 
when the drop size Is assumed constant. The reduction In the 
response peak of the spray Is due to the reduced response of the 
smaller droplets. 

Figure 9 Is the wet bulb total response for a monodlsperse 
spray as a function of co^ at 10 atmospheres pressure. The same 

conclusion can be drawn as from the non-wet bulb case. The total 
response at pressures of 1 and 100 atmospheres for both wet and 
non-wet bulb conditions give similar results. The shapes of the 
total response curves In Figs. 8 and 9 are very similar to the 
shape of the Individual response curves for A^. 

When a polydlsperse spray Is Introduced Into the combustion 

chamber, the spray Is separated Into discrete spray sets, each 

with Its own A , and the calculations proceed as before. Figure 
o 

10 Is a plot of a non-wet bulb polydlsperse spray at a pressure 
of 10 atmospheres, consisting of three spray sets with A^'s of 

1.234 X 10^®, 1 X 10^, and 1.234 x 10^. The response from each 
size group Is weighted and summed In order to obtain the total 
response. The mass fractions of the spray secs and Injection 
rates were varied, as listed In Fig. 10. The total response is 
greater than unity when the larger A^ spray sets comprise about 

30% of the injected spray mass. The results indicate, however, 
that even a very small percentage of large drops In the inject- 
ed spray can result in substantial response peaks. Findings at 
other pressures for both wet and non-wet bulb conditions are 
similar; however, for the wet bulb case, no peaks greater than 
unity are observed. 

2.6 Summary 

The response of monopropellant droplets and sprays was 
considered using an analysis based on earlier strand combustion 
work. For individual droplets, small droplets gas low or 
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Figure 8 Total Response as a Function of ui^ for a Monodisperse 

Spray at Approximately 10 atm. Pressure (Non-Wet Bulb 
Conditions) 
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Figure 9 Total Response as a Function of o)^ for a Monodlsperse Spray at 
Approximately 10 atm. Pressure (Non-Wet Bulb Conditions). 
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negative response at all frequencies. Large drops provided a 
response near unity at low frequencies, with the response de- 
clining at frequencies somewhat In excess of the characteristic 
frequency of the liquid phase thermal wave. If there was a 
temperature gradient In the liquid phase, a response peak In 
excess of unity was determined. For large droplets, a re- 
sponse peak resulting from gas phase transient effects was also 
observed, but at frequencies generally higher than the range 
normally associated with combustion Instability. 

Spray response was also considered for several specific 
cases. The response of a spray Is generally lower than that 
of the largest Injected droplet, however, even a small percent- 
age of large droplets can yield appreciable response. 

Work In this area of the Investigation has been concluded. 


3. High Pressure Spray Combustion 
3 . 1 Introduction 


Knowledge of the gasification process of Individual drop- 
lets Is vital to understanding the entire spray combustion pro- 
cess. This aspect of the problem has already been considered 
In some detail. The results have provided a capability for pre- 
dicting gasification rates and critical combustion conditions 
In environments typical of spray combustion processes. 

In order to apply the results of Individual droplet studies 
to sprays, the structure of the spray must be determined In order 
to specify the local environment of the droplets. In spite of 
Its importance, investigations of spray structure during com- 
bustion have only recently appeared due to the experimental dif- 
ficulty of this type of work. 

Chigier and coworkers (11-13) have measured the structure 
of spray combustion processes at atmospheric pressure. Gas 
temperatures and droplet sizes and velocities were determined 
within the spray tor several injector types. Onuma and 
Ogasawara (14) conducted similar measurements at atmospheric 
pressures, with additional results on species concentrations. 

These studies indicate some similarities between the struc- 
ture of spray combustion processes and turbulent gaseous dif- 
fusion flames. The similarity is most apparent when the droplets 
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are finely atomised so that droplet velocities quickly come into 
equilibrium vith local gas velocities. Under these conditions, 
models analagous to homogeneous two phase flow (negligible slip 
conditions) are suggested, although they have not been pursued 
as yet. 

At high pressures, where the density of the liquid phase 
approaches that of the gas, the condition of negligible slip 
between the two phases becomes more appropriate. If this is the 
case, the combustion of a spray at high pressures would be similar 
to that of a turbulent gas jet with greater density variations due 
to the presence of liquid in the early stages. 

Newman and Brzustowskl (15) considered the evaporation of 
carbon dioxide sprays at elevateo pressures. Experimental re- 
sults were compared with a homogeneous two phase flow model, which 
did not consider high pressure thermodynamics and employed an early 
version of the integral model for a Jet. The success of the pre- 
clctlons was limited, but it is encouraging that results were of 
the correct order of magnitude. 

Avery and Faeth (16) considered a rather different problem, 
involving the combustion of a high velocity gaseous oxidizer Jet 
submerged in a liquid metal fuel. In this case, the gas is the 
distributed phase. Since the gas has low Interla in a liquid 
surroundings, the local slip between the two phases is small and 
the homogeneous model is a good representation of the process. 

The major uncertainty involves the proper treatment of the 
enormous density variation of the flow as mixing and reaction pro- 
ceeds. The density variation was handled using a coordinate trans- 
formation and a variable density entrainment model proposed by 
Morton (17). 

The results of the procedure were very encouraging. A 
unified correlation of flame length, temperatures and velocities 
was achieved for both the two-phase system and earlier work on 
forced turbulent diffusion flames (18). Figure 11 gives an 
example of the flame length results. The gas-llquld combustion 
process Is compared with the gas-gas measurements of Hawthorne, et 
al (18) and Wohl, et al (19). Using a similar approach It was 
also possible to correlate the measurements of length of con- 
densing vapor Jets In subcooled liquids (20, 21) In the same manner 
as shown In Figure 11. 

These findings suggest that under conditions where the slip 
between the two phases Is small, the model of Ref. (16) can handle 
the large density variations of a two phase flow during a mixing 
and combustion process. The model provides information on the 
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Fig. 11 Comparison of predicted and measured flame lengths for 
gas-liquid and gas-gas turbulent jets, Ref. 16. 
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structure of the process Including temperature, density, velocity 
and species concentrations. The mean position of the diffusion 
flame Is also predicted as well as the length required for complete 
combustion. The Input variables are well defined and Include In- 
jector size. Injection momentum. Injector flow rate, ambient dens- 
ity, stoichiometry and ambient velocity. These parameters can be 
specified unambiguously for most practical combustion processes. 

The objective of the present study Is to Investigate the 
validity of the Integral-entrainment model for spray combustion 
at elevated pressures where the assumptions of the model are most 
valid. During this report period, major emphasis of the Investi- 
gation was placed on the design and fabrication of the experimental 
apparatus. 

3.2 Theory 

A complete description of the basic theoretical model Is pre- 
sented In Ref. 15, only a qualitative description will be given here. 
A sketch of the present process Is given In Fig. 12. A liquid fuel 
jet Is Injected Into a large stagnant gas region at elevated pres- 
sure. It Is assumed that buoyancy forces can be neglected. Since 
gas and liquid densities approach one another at high pressures. It 
Is assumed that fuel droplets produced by the spray are In local 
kinematic equilibrium with the gas, l.e. slip Is neglected. 

Under these assumptions, integral equations of conservation 
of mass, momentum, energy and species are written. The entrainment 
of gas by the jet Is described by Morton's variable density entrain- 
ment expression (17). The equations can be Integrated In closed 
form to yield the variation of absolute enthalpy, allowing for 
heat of formation, velocity u, and concentration. Ay, through- 
out the flow field. The problem is closed by specifying the 
equations of state for the system which relates concentrations and 
enthalpies to temperature, the presence or absence of liquid and 
density. In regions some distance outside the spray boundary, 
an ideal gas model may be used. Within the spray boundary, and 
particularly to define its extent, the high pressure thermodynamic 
models, of Refs 22-24 must be employed. 

Aside from mean quantities, the outer spray boundary and outer 
flame boundary are characteristics of interest. The interpretation 
of these factors, however, is somewhat complicated for turbulent 
Jets due to local unmixedness or property fluctuations. For ex- 
ample, the flame zone is relatively thick for turbulent flow and 
profiles of mean fuel and oxidizer concentrations overlap to a con- 
siderable extent. Therefore the outer flame limit corresponds to 
the point where all the fuel is consumed, which is beyond the 




12 Sketch of the spray combustion process 
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•qulval«nt lamliuir •tolchion«trlc flame position. Similar behavior 
also holds at the outer spray boundary. It has been found, how- 
ever, at a constant ratio allows outer boundaries to be specified 
from mean quantities for both thermodynamic and reactive processes 
(15). 


3.3 Experimental Apparatus 

Initial evaluation of the theory for the noncombusting case 
can employ existing data in the literature (15) . An apparatus has 
been designed and fabricated for tests involving the combustion 
case. A sketch of the apparatus Is presented In Fig. 13. 

The apparatus consists of a windowed pressure vessel which 
can operate at pressures In excess of 100 bar. The Internal 
volume of the chamber Is approximately 11 liters. 

A given volume of fuel Is placed In a loop In the Injector 
supply line. This sample Is forced Into the chamber by opening 
the solenoid valve between the Injector and the air supply which 
Is maintained at selected pressures higher than the test chamber. 
The spray is ignited with a hydrogen gas Jet which, is ignited by 
a hot wire. Injector passage diameters are small to reduce wall 
effects, on the order of .2-. A mm. Differential pressure gages 
are used on the fuel feed and hydrogen supply, an absolute pres- 
sure gage Is used for the chamber pressure. 

The process Is photographed with a motion picture camera, 
using a dark field for flame positions, or back-lighting and 
filters for spray positions. Measurements can also be made of 
gas temperatures in the Jet using thermocouples recorded by an 
oscillograph. 

The time sequencing of injection, ignition, shutdovm, etc., 
is controlled automatically using a programmable timer. Initial 
work is considering methanol and heptane as fuels, due to avail- 
ability of developed thermodynamic models for the high pressure 
characteristics of these materials (22-24) . 

3.4 Summarv 


A theory for high pressure spray combustion has been obtained 
by modifying an analysis developed for gas-liquid and gas-gas 
combustion and condensation processes. Modification of earlier 
high pressure thermodynamic analysis is in progress for use in pre- 
dicting spray boundaries with this model. Initial model evaluation 
is being undertaken using existing data in the literature. 
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Fig. 13 Sketch of the high pressure spray combustion apparat 
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An apparatus has baen designed and fabricated to allow ob- 
servation of combusting fuel sprays at high pressures. The 
apparatus provides measurements of spray shadowgraphs, dark 
field flame photographs and gas temperatures in the region 
downstream of the spray. Assembly of the apparatus is nearing 
completion and testing will begin early in the next report 
period. 
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